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ABSTRACT 

We study the evolution of the Star Formation Rate Function (SFRF) of massive 
(M^ > 1O^°M0) galaxies over the 0.4 < z < 1.8 redshift range and its implications 
for our understanding of the physical processes responsible for galaxy evolution. We 
use multiwavelength observations included in the GOODS-MUSIC catalogue, which 
provides a suitable coverage of the spectral region from 0.3 to 24 fim and either spectro- 
scopic or photometric redshifts for each object. Individual SFRs have been obtained 
by combining UV and 24 /ini observations, when the latter were available. For all 
other sources an "SED fitting" SFR estimate has been considered. We then define 
a stellar mass limited sample, complete in the M^, > lO^^M© range and determine 
the SFRF using the l/Vmax algorithm. We thus define simulated galaxy catalogues 
based on the predictions of three different state-of-the-art semi-analytical models of 
galaxy formation and evolution, and compare them with the observed SFRF. We show 
that the theoretical SFRFs are well described by a double power law functional form 
and its redshift evolution is approximated with high accuracy by a pure evolution of 
the typical SFR (SFR*). We find good agreement between model predictions and the 
high-SFR end of the SFRF, when the observational errors on the SFR are taken into 
account. However, the observational SFRF is characterised by a double peaked struc- 
ture, which is absent in its theoretical counterparts. At z > 1.0 the observed SFRF 
shows a relevant density evolution, which is not reproduced by SAMs, due to the well 
known overprediction of intermediate mass galaxies at z 2. Semi- analytical models 
are thus able to reproduce the most intense SFR events observed in the GOODS- 
MUSIC sample and their redshift distribution. At the same time, the agreement at 
the low-SFR end is poor: all models overpredict the space density of SFR ^ IMq/ut 
and no model reproduces the double peaked shape of the observational SFRF. If con- 
firmed by deeper IR observations, this discrepancy will provide a key constraint on 
theoretical modelling of star formation and stellar feedback. 

Key words: galaxies: evolution - galaxies: fundamental parameters - cosmology: 
observations 



1 INTRODUCTION 

The evolution of the star formation rate (SFR) over the 
cosmic time is a fundamental constraint fo r every theory 
of galaxy form ation and evohition (see e.g. [Hopkins' 2004; 
iHopkins &: Beacom .2006 ). Estimating SFRs for individual 



galajcies is a complex task, owing to the uncertainties in- 
volved in the reconstruction of this quantity from obser- 
vational data. It is widely accepted that dusty molecular 
clouds are the main sites for star formation: this implies 
that newly born stars are subject to significant dust atten- 
uation, until they are able to escape or disrupt their parent 
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cloud. Young OB stars emit a considerable amount of energy 
in the restframe ultraviolet (UV) band, which has thus been 
considered a key waveband for recovering the instantaneous 
SFR in external galaxies. However, dust absorbs UV pho- 
tons, heats up and re-emits this energy as thermal radiation 
in the Infrared (IR) bands. For this reason, SFR estimates 
based on UV luminosity include correction facto rs to ac- 
coun t for dust attenuation and re-emis sion (Calzetti et al.l 
11994 iKennicutd Il998l : iBell et al.l I2003D . Dust attenuation 
is particularly relevant for star forming galaxies, where the 
dust emission peak is the dominant compone nt of the galac- 
tic sp ectral energy distribution (SED, see e.g. ICalzetti et ahl 
|2000| ). n ot to mention the extr eme cases of sub- mm galaxies 
(see e.g. IChataman et al.ll2004l ). The study of the cosmic IR 
background shows that the global energy emitted by galaxies 
in the IR is comparable to the direct starlight emis sion, de- 
tectable in the optical (see e.g. iLagache et al]|2005l '). clearly 
showing that a significant fraction of star formation activity 
is expected to be heavily extinguished and detectable only 
in the IR. 

These uncertainties become more relevant, when SFR 
estimates for galaxies covering a wide redshift range are 
considered, since both galaxy physical properties and 
dust properties are expected to evolve w i th cosmic time 
("Maiolino et al."2004'; 'Fontanot et al."2009a|; [ Gallerani et al.l 
[2010. : .Fontanot fc Somerville 2011 ). Combining UV informa- 
tion with supplementary information from direct observa- 
tions in the IR region has thus be en proposed as the best 
tool t o account for the total SFR (Ke nnicutti 19981 : iBell et al] 
In particular, observations at 24fim have shown to be 
extremely useful for esti mating the global IR luminosity and 
the i nstantaneous SFR l|Papovich et al.ll2007l : [Santini et al.l 
I2OO9I ). Thanks to the advent of the Herschel Space Obser- 
vatory, we will be able to constrain total IR luminosities, by 
directly sampling the peak of the thermal emission. In fact, 
it has been recently shown bv lRodighiero et al.l (|2010l ). that 
the combination of 24^m data with Herschel observations 
at longer wavelengths represents a very promising tool for a 
more accurate determination of SFRs of individual galaxies 
up to 2 3. 

Nowadays, our view of the evolution of galaxy prop- 
erties has substantially changed thanks to the advent of 
multiwavelength surveys. If the spectral sampling is fine 
enough, these catalogues can be used to infer the physical 
properties also in the absence of spect roscopic information . 
In recent years sev eral groups (see e.g. iFontana et al.|[2004l : 
IPanter et allbOOTl ') developed a number of "SED fitting" al- 
gorithms: these codes compare the available photometry for 
individual sources with synthetic SED libraries and the best- 
fit model template is chosen by means of a minimisation. 
The resulting estimate for redshift (the so-called photomet- 
ric redshift), stellar mass (M^), dust extinction and SFR are 
the most widely used results of this procedure and have been 
of fundamental importance for our understanding of galaxy 
evolution. 

The common interpretation of these results is connected 
to the so-called "downsizing" scenario, in which the star for- 
mation shifts from high mass t o low mass galaxies as redshift 
decreases (first introduced bv lCowie et al. 199t ) . This pic- 
ture has been recently revised bv lFontanot en iT ll2009bl ) by 
showing that the typical errors associated with the estimate 
of the physical quantities has to be taken into account, when 



comparing the observed evolutionary trend s with the pre- 
dictions of theoretical models. In particular, iFontanot et ahl 
(2009h) concluded that the discrepancies seen between mod- 
els and data for massive galaxies are not significant, if model 
predictions are convolved with typical observational errors, 
while the strongest discrepancies between model and data 
are seen for low-to-intermediate mass galaxies. 

The importance of the study of the redshift evolution of 
the cosmic SFR for galaxies at a given mass r ange has been 
widely recognised by a rmmber of authors (iNoeske et al.l 
2007; Elbaz et al. 2007; Zhen g et al. 2Qp7; I)rorv fc Alvar^ 
20081 : IDunne et al.ll2009i : ISantini et all 120091 : icilbank et all 
2011f). These results and the evolution of the stellar mass 



density provide fundamental information about the evolu- 
tion of the global process of galaxy formation. The latter 
approach is best complemented by the analysis of the stel- 
lar mass function (i.e. the volume density of galax;ies as a 
function of stellar mass), which allows us to characterise 
galaxy evolution as a function of both redshift and stel- 
lar mass. Despite the wealth of information already avail- 
able in the literature about the stellar mass function and 
its redshift evolution, very little is known about the corre- 
sponding Star Fo rmation Rate Fu nction (SERF hereafter) 
and its evolution: iBell et all |20o3) studied its evolution in 
the 0.2 < z < 1.0 redshift r ange in the COMBO-17 sur- 
vey, [Schimingvic^^eF|al] (|2007l ) presented the di s tribut ion of 
SFRs in the SDSS survey and iBothwell et"al] l|201ll ) con- 
sidered the distribution of SFRs in the Local Volume (de- 
fined as an 11 Mpc radius around the Milky Way) combining 
GALEX and Spitzer observations. This is mainly due to the 
uncertainties involved in the SFR estimate, but also to the 
difficulties in correctly accounting for the completeness cor- 
rection needed, especially at the low-SFR end. 

In this paper, we focus on a new treatment of the 
GOODS-MUSIC data, by focusing on a mass complete sam- 
ple of galaxies. In this way, we are able to fully characterise 
the various levels of star formation in these objects and fol- 
low their redshift evolution. Moreover, since the strongest 
starbursts are hosted by galaxies above the chosen mass 
limit, the high-SFR end of the resulting SERF provides fun- 
damental information about the redshift evolution of this 
particular galaxy class. Therefore, this paper represents a 
complementary study of the evolutio n of SFRs in galaxie s 
of different stellar mass with respect to lSantini et al. (|2009h . 
to which we refer the reader for the analysis of the evolution 
of the cosmic SFR and specific SFR. A detailed comparison 
of the same catalogue with other estimates for the evolution 
of SFR as a function of red shift and stellar mass has been 
presented in lFontanot et al] (f2009b,) (see references herein). 
In the present paper we focus on the determination of the 
SFRF and compare results with the predictions of state- 
of-the-art semi-analytical models of galaxy formation and 
evolution. 

This paper is organised as follows. We present the 
GOODS-MUSIC catalogue in sec. [2] and the theoretical pre- 
dictions in sec. [S] we present our results in sec. [4] finally we 
discuss our conclusions in sec. [S] In the original GOODS- 
MUSIC papers a Salpeter initial mass function (IMF) is 
mostly used; throughout this work, we assume a Chabrier 
IMF and we rescale all stellar mass and SFR estimates, ac- 
cordingly. We also adopt the ACDM concordance cosmolog- 
ical model {Ho = 70km/ s /Mpc, Qm = 0.3 and Q.x = 0.7). 
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2 DATASET 

Our dataset is ba sed on the updated version of GQQD - 
MUSIC catalogue (|Grazian et al.ll2006l : ISantini at al.ll2009^ . 
GOODS-MUSIC is a multiwavelength catalogue, with pho- 
tometry in 15 bands, from 0.35 to 24 /im, covering an area 
of ~ U3.2arcmin^ of the GOODS-South field and listing ~ 
15000 sources. Main source detection has been performed 
on the z image, but faint z-band objects detected either 
on the K image or at 4.5 fim have also been included in 
the catalogue. The GOODS-MUSIC catalogue has been also 
cross-correlated with available spectroscopic catalogues, and 
a spectroscopic redshift has been assigned to ~ 12% of 
sources. For each source in the catalogue, its physical param- 
eters, such as Mi,, SFR and the photometric redshift, have 
been estimated thr ough a "SED fitting" algorithm (see e.g. 
ISantini et ahllioogl and reference herein) ; for the purposes of 
minimisation, the available broad band photometry up 
to rest-frame A < 5.5/im has been consi dered. The reference 
synth etic SED library is defined using iBruzual fc CharlotI 
l|2003l ) synthetic models and assuming exponentially declin- 
ing star formation histories. The relevant parameters for 
the SED synthesis (such as dust content, metalicity, age, 
timescale for star formation histories, etc) are defined on 
a grid in order to cover a wide range of possible values 
and combinations. The discreteness of SED generation im- 
plies that the final estimates for the physic al parameters are 
subje ct to relevant degeneracies (see e.g. iMarchesini et all 
I2OO9I for a discussion focused on the effect of changing the 
physical ingredients in the adopted stellar population mod- 
els). In particular, theoretical models predict star formation 
histories for individual galaxies, which grossly differ from 
any assumed smooth analytical s hape: this has impo rtant 
implications for the M* estimate l|Stringer et al.ll2009l ). By 
comparing the available spectroscopic redshifts (zspcc) with 
the corr esponding estimates from the SED fitting algorithm 
(zphot), ISantini et al.l (|2009l ) computed the average abso- 
lute scatter between spectroscopic and photometric redshift 
They found a{Az) = 0.06 for the whole 
= 0.043 when only the brightest galaxies 
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sample and a{Az) 
{z < 23.5) are considered. 

A parallel estimate of the star formation rate is pro- 
vided by the combination of IR and UV luminosities. For 
sources in the GOODS-MUSIC sample, [Saiitini et al.l (|2009l ') 
combined the rest-frame UV luminosity derived from the 
SED fitting and the integrated luminosity between 8 and 
1000 micron, obtained by fitting 24 micron fluxes with 
Dale fc Helod (20021) t emplates. The correction provided by 
Papovich et all (|2007l ) has been applied on the resulting 



SFRs in order to account for the overestimation of the to- 
tal IR luminosity at high redshifts/IR luminosities when 
extrapolating from 24 fim fluxes (e.g. ^Elbaz ct al. 2010). 
ISantini et al.1 (120091 ) explicitly compared these SFR esti- 
mates with the SED fitting results, for sources with both 
estimates available. They found an overall good agreement, 
but also a systematic trend at all redshift, with IR-based 
estimates exceeding the fit-based ones as the star forma- 
tion mteJri^reaseSjaBdv of SFR (see 
also lPapovich et alllioOTl : iNordon et al.ll2010l ). In the follow- 
ing of this paper, unless explicitly stated, for each object in 
the sample we consider IR-based SFRs when available, and 
"SED fitting" estimates elsewhere. 



In order to provide a reliable estimate of statistical 
quantities like the stellar mass function and the SFR func- 
tion, a careful analysis of observational biases and selection 
effects is of fundamental importance for a complete census 
of all galaxies. Those corrections are critical for the SERF, 
given the complex combination of multiwavelength informa- 
tion needed, including both the selection of galaxies in the 
AT-band, fundamental for the definition of a mass-limited 
sample, and the SFR estimate from 24/im for direct de- 
tection (if available) or UV and optical data for the SED 
fitting procedure. In order to reduce the uncertainties in- 
volved in our analysis, we decide to focus on a subsample of 
the GOOD-MUSIC catalogue, defined by considering only 
sources with a stellar mass estimate M* > 10^^ Mq. The 
GOODS-MUSIC catalogue is expected to be highly com- 
plete in this mass range over the whole 0.4 < z < 1.8 red- 
shift interv al, corresponding to _ft "-band magnitudes brighter 
than ~23 (|Fontana et al.l |2006| ). Moreover, for these ob- 
jects, either 24/im fluxes or detailed multiwavelength pho- 
tometry are available, and provide reliable estimates for 
SFR> O.OIMq /yr (Santini et al. 2009). Objects with docu- 
mented AGN activity (i.e. spectroscopically determined or 
X-ray sources) have been removed from the final sample. 
Compton thi ck candidates, sel ected on the basis of their very 
red spectra (|Fiore et al.ll2008l ). have been removed as well. 
In the redshift range of our interest there are only 4 removed 
objects with 1.4 < z < 1.8 (0.5% of the considered sample). 
We test that our conclusions are robust against the inclusion 
of these objects. Our final sample consists of 740 sources. 

It is also worth mentioning the known overdensities 
at z 0.7 and z ~ 1.6 in the GOO DS fields, as well 
as the known underdensity at z ^ 0.9 (|Gilli et al.l l2003l : 



ICastellano et aL l l2007l:ISalimbeni et al.ll2009D . On going sur- 
yeys like COSMOS Is^viUe et al.1 12OO7I) and CANDELS 
(|Grogin et al.ll201ll : Koekemoer et al.ll201ll ). wiU be able to 
reduce the effect of cosmic variance on the SERF determi- 
nation, and study its environmental dependencies. 



3 MODELS 

In order to understand the physical implications of our re- 
sults, we consider the prediction of the recent implementa- 
tio ns of 3 state - of-the -art semi-an alytical models, namely , 
the IWang et al.l ||2008|) model, the ISomerviUe et al.1 (I2008D 
and MORGANA (|Monaco et al.l |2007| ). The predictions of 
these 3 models ha ve been extensively compared in a number 
of re cent papers l|Kimm et al.ll2009l : iFontanot et al 1 l2009bl , 
I2OIII ) and we refer the reader to these publications and the 
original papers for all the details of the galaxy formation 
and evolution modelling. In these models, the evolution of 
the baryonic component is followed by means of an approx- 
imate description of the physical processes at play (such 
as gas cooling, star formation, stellar feedback. Black Hole 
growth and AGN feedback) and of their interplay with grav- 
itational processes (i.e. dynamical friction, tidal stripping 
and two-body mergers) , linked to the assembly of the large- 
scale structure of the Universe. These "recipes" include a 
number of parameters which are usually fixed by comparing 
model predictions with a set of low-redshift observations. 
Despite their simplified approach, SAMs have turned into a 
flexible and widely used tool to explore a broad range of spe- 
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Figure 1. The redshift evolution of tiie SFRF in semi-analytical 
models. Crosses refer to the intrinsic predictions of the models, 
while solid lines show the analytic fit (pure SFR* evolution) of the 
theoretical SFRF. In each panel, dashed line reports the position 
of the 2 = SFRF. 



cific physical assumptions, as well as the interplay between 
different physical processes. All models we consider resolve 
galaxies with M* > 10^ Mq, and they reproduce fairly well 
a number of observational constraints on the evolution of 
the SFR, both as a global quantity (i.e. the cosmic star for- 
mation rate), and for individual galactic population (i.e. the 
fraction of passiv e galaxies in the local Universe, see e.g. 
iKimm et allboogl '). 

Also for the theoretical predictions, we consider purely 
stellar mass selected samples (M* > 10^°Mq) and, as for the 
GOOD-MUSIC catalogue, we excluded from the analysis all 
objects with SFR< O.OlMQ/i/r. After the sele ction proce- 
dure w e are thus left with 385 187 galaxies for th e Wang et al.l 
l|2008l ) model, 342795 for the lSomerviUe et al.l (|2008l ) model 
and 249823 galaxies for MORGANA. In fig. [T] we show the 
binned SFRFs for each model (crosses). In order to quantita- 
tively estimate the redshift evolution of the predicted SFRF 
we apply to mode l predi ctions a fitting procedure similar 
to iFontanot et al.l l|2007ar ). We determined best fit analyt- 
ical solutions to the theoretical predictions assuming both 
a Schechter and a double power law shape for the SFRF. 
We use the corresponding number densities to Monte Carlo 
generate independent data points in the z-SFR space and we 
calibrate the best fit parameters for the SFRF and its evolu- 
tion by means of a binning procedure. Our results show 
that the double power law shape systematically scored bet- 
ter results than the Schechter function in terms of x'^ (with 
typical values < 1.0). This is mainly due to the extra 
degree of freedom in the high-SFR-end slope, with provides 
a better description of this region with respect to the expo- 
nential cut-off of the Schechter function, even if it is worth 
stressing that the Schechter shape is not formally ruled out 
by our analysis. Anyhow, in the following we thus assume a 
SFRF in the form of a double power-law: 



${SFR) = 



^*{SFR*) 



(1) 



(SFR/SFR*)-" + {SFR/SFR*)-i^ 

where SFR*, a, P are free parameters. In order to model 
the redshift evolution of the SFRF we fix the values of the 
4 parameters at their 2 = values and we consider different 
combinations for their redshift evolution. We consider both 
a pure evolution in SFR* (PSE) of the form: 

SFR* = SFRUoC^ + ■z)'"'"''' (2) 
and a pure density evolution (PDF) of the form: 

$* = <l>:^oe**^'+"' (3) 

where fcsFR and fc$ represents the evolutionary parameters. 
Our results point out that the redshift evolution of the the- 
oretical SFRF is in general well described by a pure evo- 
luti on in SFR*. In fig . [TJ we show the best fit results for 
the IWang et al.l (|2008l ) model as a solid red line, while the 
dashed magenta line represents the corres ponding z = 0.1 
best fi t SFRF. As for morgana and the ISomerville et al.1 
l|2008l) model, the fit at the low-SFR end is improved (by 
a factor of 2 in terms of x^) if the pure SFR* evolution 
(which is still statistically acceptable) is combined with an 
evolution of the slope a of the form: 



a = ao + kaZ 



(4) 



where ao = a{z = 0) and ka is the evolutionary parameter. 
We conclude that a double power law is a good represen- 
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Table 1. Best fit parameters for the evolution of the SFR function as predicted by semi-analytical models. 



SAM 


**(xl0~2) 


SFR* 


ksFn a 


/9 


ao 






Wang et al. (2008) 
Somerville et al. (2008) 

MORGANA 


0.47 ±0.04 
0.46 ± 0.05 
0.41 ±0.03 


0.29 ± 0.14 
0.39 ±0.12 
0.70 ±0.17 


2.10 ±0.64 0.84 ±0.10 
1.99 ±0.72 — 
1.23 ±0.42 — 


-2.03 ±0.27 
-1.84 ±0.24 
-1.95 ±0.22 


0.29 ± 0.07 
0.26 ± 0.03 


0.23 ±0.11 
0.15 ±0.04 


0.32 
0.44 
0.29 



tation of the SFRF as predicted by the theoretical models 
at J3 < 2 and for SFR> O.OlMg/yr. We present the best fit 
parameters corresponding to the above equations in tab. [U 
It is worth stressing that we obtain very low value even 
if our fits tend to break down at the highest redshift and 
lowest SFRs. This is due to the high number of object in 
the model samples, which thus provide a very accurate fit 
of the knee of the SFRF and its evolution. 

In order to compare the theoretical predictions with the 
results from the GOODS-MUSIC catalogue, we repeat the 
fitting procedure after convolving the theoretical predictions 
with an estimate of the systematic error on the SFR determi- 
nation, to fully account for the typical observational error in 
the catalogue ( Santini etjd. 2009). We assume a log-normal 
error dis tribution for the SFR with an amplitude of 0.3 dex 
(see also lFontanot et al.ll2009bl ). We repeat the above analy- 
sis and we determine that a pure SFR* evolution of a double 
power law is still a good description (i.e. consistent val- 
ues) of the overall redshift evolution of the predicted SFRFs. 
The most noticeable difference in the best fit parameters is 
the systematic increase of the (5 parameter by as much as 
20%. 



4 RESULTS & DISCUSSION 

We first compute the binned SFRF for massive galaxies by 
means of the standard l/Knax formalism, and we show the 
resulting SFRF in fig. [2] as open circles (cyan circles mark 
the bin with less than 5 sources). In order to take into ac- 
count the observational uncertainties, we repeat the anal- 
ysis 10000 times after convolving every individual SFR es- 
tima te with the corresp onding error on SFR and redshift 
(see ISantini et al.l |2009| for more details on the error es- 
timate). Errorbars in fig. [2] refer to the 1-a variance over 
these 10000 realisations. We collect our final estimate of the 
SFRF and associated errors in tab[2l In the fig. [2] we also 
include the previous esti mate for the SF RF at 2 > from 
the COMBO-17 survey (jBell et all boOTl 'l . It is worth not- 
ing that the COMBO-17 results are obtained by considering 
all galaxies in their sample, with no cut in stellar mass. 
Our estimates for the SFRF are in good agreement with 
thei_BcU et al. (2007) result at the high-SFR end, where we 
expect the massive galaxies to dominate the space density. 

Fig. [2] clearly show that the COMBO-17 data are not 
deep enough to explore the low-SFR end with the same ex- 
tent as in the GOOD-MUSIC catalogue, moreover the low- 
SFR end of the SFRF shows a very peculiar feature, i.e. 
a double peaked shape. It is really tempting to interpret 
this feature in terms of the well known bimodality in the 
colour (SFR)-mass diagram, and/or connect those gala:xies 
with the intermediate sources populating the so-called green 
valley. It is also worth noting that a populations of rela- 



Table 2. SFR Function in the GOODS-MUSIC Catalogue. 



log{SFR) Log{<S>{SFR)) 



0.4 < 2 < 0.6 



-1.75 ± 0.25 
-1.25 ± 0.25 
-0.75 ±0.25 
-0.25 ± 0.25 
+0.25 ± 0.25 
+0.75 ± 0.25 
+1.25 ± 0.25 
+ 1.75 ± 0.25 



0.6 < 2 < 0.8 



-1.75 ± 0.25 
-1.25 ± 0.25 
-0.75 ± 0.25 
-0.25 ± 0.25 
+0.25 ± 0.25 
+0.75 ± 0.25 
+ 1.25 ± 0.25 
+ 1.75 ± 0.25 



0.8 < 2 < 1.0 



-1.75 ± 0.25 
-1.25 ± 0.25 
-0.75 ± 0.25 
-0.25 ± 0.25 
+0.25 ± 0.25 
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Figure 2. The evolution of the SFRF for M > IO^^A/q galaxies in the GOOD-MUSIC sample. In each panel, the horizontal cyan long- 
dashed line represents the space density corresponding to one object in the GOODS-MUSIC volume, while the yellow shading highlights 
the the SFR range covered by 2i^ra observations in each redshift bin. Open circles refer to our SFRF estimate from the GOODS-MUSIC 
survey: in each subpanel, cyan solid circles mark the bins with less than 5 galaxies; errorbars represent the 1-cr confidence levels (see text 
for more details). Light grey filled circles refer to the SFRF derive d only from the results of the SED fitting algorithm. Purple diamonds 
refer to the SFRF estim a te from the COMBO- 17 llBell et al.ll200'1) survey. In all panels, dashed red, dot-dashed green and solid blue lines 
refer to the lWane et al.l ||2008| ) . [Somerville et al.] i 20081) and MORGANA predictions, respectively. SAM predictions are convolved with an 



Gaussian error on the SFR (width 0.3 dex, see text). 



tively massive galaxies with very low levels of star formation 
has been already reported at lower redshifts (|Wilman et al.l 
l2008l : IScliawinski et aLllioioh . 

However, a word of caution is necessary. We have lim- 
ited information about individual 24pm fluxes for galax- 
ies in our mass-selected sample, and typically only for 
SFR> IMQ/yr. This implies that the secondary peak at 
SFR< IM^/yr is populated by galaxies whose SFR esti- 
mate comes from the SED fitting algorithm. In order to test 
the robustness of our result, we repeat our analysis consid- 
ering only SFR estimates coming from the SED fitting. We 
find no significant change in the shape of the resulting SFRF 
(light grey filled circles in fig. [2]): in particular, the double 



peaked behaviour persists. Our main conclusions are there- 
fore insensitive to t he ori gin of SFR estimate, in agreement 
with ISantini et all (|2009l ') results (i.e. the 24^im-based SFR 
estimates are in reasonable agreement with SED fitting re- 
sults for sources with both estimates available). 

Nonetheless, we cannot completely exclude the hypoth- 
esis that the gap between the two peaks is an artifact of 
the SED fitting procedure (i.e. piling up of low-SFR galax- 
ies in the secondary peak, due to the coarse grid-sampled 
properties of the synthetic SED library). We plan to study 
the efficiency of the SED fitting algorithm in the estimate of 
SFRs in the range O.IMq /yr < SFR < WMq /yr as a part of 
a larger programme, comparing the global performances of 



© 0000 RAS, MNRAS 000, 000-000 



SFR function from GOODS-MUSIC 7 



our algorithm in recovering the physical properties of model 
galaxies with complex star formation histories, starting from 
synthetic photometry extracted from SAM simulated cata- 
logues. A simila r procedure has been recently employed by 
iLee et alj (|2009l ) to a sample of high-redshift model galax- 
ies: they found that SED fitting underestimates SFR in their 
sample and attributed this behaviour to to the different star 
formation histories predicted by SAMs with respect to those 
assumed in the reference synthetic SED libraries. 

As a preliminary test we consider a small sample of 
5000 model galaxies with 0.4 < z < 1.8 taken from the same 
MORGANA realisation and uniformly distributed in SFR. For 
the purpose of this test we do not limit the sample to 
M* > W^'^Mq model galaxies. For each object, we compute 
a synthetic SED using the Radiative- Transfer code GRASIL 
ISilva et al.lll998l . see also [Fontanot ct al. 200713). We then 
apply the SED fitting algorithm to the corresponding syn- 
thetic photometry. In fig. [3] (upper panel) we compare the 
original, flat, distribution of SFR in our sample (black his- 
togram) with the distribution of reconstructed SFRs (blue 
shaded histograms). This test suggests that some depletion 
effects in SFR distribution may be effective at SFR below 
our fiducial limit (SFR O.OlMo/yr); however, the SFR 
distribution in the region corresponding to the secondary 
peak and the valley (roughly identified by the green arrows) 
does not seem to be affected by any obvious systematic ef- 
fect. We also repeat the same analysis on a subsample of 
Mi, > model galaxies, (lower panel), whose distri- 

bution is more similar to the theoretical SERF. Again we 
obtain no obvious systematic effect in SED fitting recon- 
struction. We also check that similar results hold if we split 
the sample in smaller redshift bins, but with less statistical 
significance, due to the reduced dimensions of the subsam- 
ples. 

In fig. [2] we also show the predictions of our three SAMs 
as coloured lines: the level of agreement between the 3 mod- 
els is striking, given the very different star formation and 
stellar feedback schemes adopted. In particular we stress 
that the predicted redshift evolution of the typical SFR* 
is remarkably similar between all models. This is due to the 
external constraints (i.e. the evolution of the cosmic SFR) 
used to calibrate the models. The strongest discrepancies are 
seen at the low-SFR end of the SERF and clearly show the 
importance of using constraints based on the SERF, in or- 
der to disentangle between different theoretical approaches 
to star formation and stellar feedback. 

The comparison between theoretical predictions and 
the observational results from the GOODS-MUSIC survey 
shows a remarkably good agreement at the high- SFR end, 
where the available 24/im observations are deep enough to 
constrain its shape and evolution out to z ~ 1.8 (yellow 
shaded region in fig. [2|. Therefore we conclude that SAMs 
provide a satisfactory description for the space density of 
galaxies with SFR> lOAlQ/yr at moderate redshifts. It 
worth stressing that accounting for observational errors is 
essential for the correct interpretation of the comparison be- 
tween model and data: the error-free intrinsic SEREs would 
underpredict the space density of the strongest starbursts 
by a factor of a few. 

Observations at 24^m are also able to constrain the 
position of the peak of the SERF out to z ~ 1 (and give 
some indications for its position at higher redshifts). In or- 
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Figure 3. Test on the SED fitting algorithm. Upper panel: com- 
parison of the SFR distribution in a sample of 0.4 < z < 1.8 
MORGANA model galaxies (solid black histogram), with the corre- 
sponding distribution of reconstructed SFRs, obtained applying 
the SED fitting algorithm to synthetic photometry (blue hatched 
histogram ). Red l ine re presents the limiting SFR considered in 
the Santin i et al paper, while green arrow bracket the re- 

gion containing the main features of the observed SERF. Lower 
Panel: same as upper panel, but for a sample of M* > W^^Mq 
model galaxies. 

der to get more insight on the evolution of the typical SFR 
of galaxies at different redshifts, we repeat on the data the 
same fitting procedure described in sec. [S] Given the double 
peaked shape of the observed SERE both a double power 
law and a Schechter function are no longer a good represen- 
tation of the SERF over the whole SFR range. Therefore, we 
chose to repeat the fitting procedure both on the reference 
sample {SFR > O.OIAIq /yr), and on a smaller sample con- 
taining only sources with SFR > IMq/dt (i.e. focusing on 
the primary SERF peak). This choice has no major impact 
on our conclusions; as expected only the low-SFR-end slope 
a show a relevant variation between the two samples. We 
collect our final results in tab. [3] The results show that the 
SERE, as depicted by the GOODS-MUSIC survey, shows a 
remarkable density evolution (PDE) in the 0.8 < z < 1.8 
range, and a milder SFR* evolution (PSE) over the whole 
redshift range. The theoretical SERE is able to reproduce 
the SFR* evolution, although it peaks systematically at 
lower typical SFR than the observational constraints, but 
its density evolution is negligible. This represents a first rel- 
evant discrepancy between models and data and it is due 
to the well documented overprediction of the space density 
of 10^°Mq < < 10^^ Mq gala xies at 1.0 < z < 1.8 
seen in the models and in the data : iFontanot et al.l l|2009bD 
showed that SAMs predict for these galaxies space densities 
very similar to the local value out to z ~ 1.8, while their 
observed stellar mass function show a remarkable evolution 
(their fig.l). This implies that the negligible density evolu- 
tion predicted for the SERE at z > 1.0 is related to the too 
efficient formation of interm ediate mass galaxies in SAMs 
(see also lLo Faro et al.l|2009l ). 

The most interesting discrepancies between models and 
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Figure 4. Upper panel: Specific SFR versus Af* Lower panel: Rest-Frame colour-magnitude diagrams. In each panel the yellow dots 
refer to the _fC-limited GOODS-MUSIC sample; red dots mark sources with M* > 101°Mq and SFR > O.OliWe/yr; blue crosses refer 
to galaxies belonging the secondary peak on the SFRF at each redshift bin, while green circles refer to galaxies with SFR intermediate 
between the two peaks. In the lower panel, dark red and blue lines refer to the mean U — V colours alon g the red sequence and the blue 
cloud, while the dark green solid line reproduces the divider between the blue and red populations (see ISalimbeni et al.ll20o3 for more 
details) . 
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Table 3. Best fit parameters for the evolution of tiie SFR function in the GOODS-MUSIC Catalogue. 



Model 


4>*(xl0~2) 


SFR* 


fcsFR 




a 


/3 










Log (SFR) > 


-2.0 








PDE 
PSE 
PDE+PSE 


0.37 ±0.11 
0.74 ±0.20 
0.88 ±0.16 


1.13 ± 0.13 
1.30 ±0.07 
0.81 ± 0.11 


0.70 ±0.58 
2.07 ±0.44 


-1.09 ±0.60 
-1.22 ±0.22 


unconstr. 
unconstr. 
unconstr. 


-1.77 ±0.30 
-1.85 ± 0.29 
-1.43 ± 0.31 


2.7 
1.6 
1.2 


Log{SFR) > 0.0 


PDE 
PSE 
PDE+PSE 


0.30 ±0.05 
0.58 ±0.16 
0.55 ±0.05 


1.02 ±0.19 
1.31 ±0.11 
0.54 ±0.15 


0.93 ±0.35 
2.09 ±0.38 


-0.72 ±0.23 
-0.56 ±0.11 


unconstr. 
unconstr. 
unconstr. 


-1.84 ±0.21 
-1.48 ± 0.24 
-1.43 ± 0.24 


1.5 
0.9 

0.6 



data are indeed seen at the low-SFR end: all SAMs predict a 
featureless low-SFR end, whose slope (and its redshift evo- 
lution) depends on the details of the treatment of SFR and 
stellar feedback. The most relevant difference between data 
and model prediction is the fact that all models overpredict 
the abundance of SFR~ IM^/yr sources. It is also interest- 
ing to notice that no model reproduces the double peaked 
feature of the observed SFRF. 

There are at least two possible interpretations of our re- 
sults. First (scenari o A), as it has been suRge st ed by many 
auth ors (see e.g., iDe Lucia fc BlaizotI l2007l : IWang et al.l 
l2007l 'l. SAMs predict too low SFRs for most of the galax- 
ies in the mass range of our interest: they have SFRs well 
below the detection limit of our sample, an d, overall, they 
are to o passive at intermediate redshifts (|Fontanot et al.l 
l2009bl ). The three models indeed predict a secondary peak 
in their SFR distribution, and it consists of completely pas- 
sive g alaxies (SFR= 0, see e.g. figure 4 in iFontanot et al.l 
l2009bl ). In this scenario, we could interpret the observed sec- 
ondary peak (at low, but still detectable SFRs), if confirmed 
by deeper observations, as an evidence that models are not 
able to reproduce the correct space density for galaxies mi- 
grating from the blue cloud to the red sequence (see also 
ICassata et al, 2007), and we could interpret it as an indica- 
tion that galaxies in the GOODS-MUSIC sample move in 
the colour-magnitude diagram at a slower pace with respect 
to theoretical expectations. This scenar i o is al so consistent 
with the recent findings of lLagos et al.l (|201ll ): they revis- 
ited the definition of the star formation law in the context o f 
the GALFORM SAM (|Baugh et al.ll2005l : ISo^er et al.l 120061 ') 
taking into account recent results on the relation between 
star formation a nd the atomic and molecular gas content of 
galaxies (see e.g. iBlitz fc Rosolow skv 20061: iKrumholz et al.l 
[2009). They showed that the properties of galaxies with low 
SFRs (and their redshift evolution on the M* versus SFR 
plane) are particularly sensitive to the details of the assumed 
star formation law, and suitable modifications to the stan- 
dard prescriptions may indeed lead to bimodal SFR distri- 
butions in the direction needed to explain our results. 

An alternative scenario (scenario B) requires small 
SFR retriggerin g events in already red galaxies (see e.g. 
iTreu et ahl 20051 ) : these events might be linked to the dy- 
namical evolution of galaxies in dense environments (merg- 
ers, interactions). Therefore this also implies a dependence 
of these secondary SFR events on the environment, which 
in principle could b e accessible by larger surveys (see e.g. 
ICucciati et al.|[2010l and references herein). Another possi- 



ble origin for small and late star formation episodes in mas- 
sive galaxies might be related to the duty cycles associate d 
to the "radio- mode" AGN activity (|Fontanot et all [2OI1I ') : 
if the injection of energy from the central Black Hole into 
the surrounding Dark Matter halo is not a continuous event, 
some cold gas might still be allowed to infall onto massive 
galaxies. 

In order to get more insight in the properties of galaxies 
responsible for this feature, we consider the distribution of 
the GOODS-MUSIC galaxies both in the specific SFR versus 
Mi, and in rest-frame colour-magnitude diagram (as recon- 
structed from the SED fitting algorithm) in fig |4] In both 
panels, we mark the M* > lO^^M© galaxies as red squares; 
blue crosses and green circles refer to the galaxies belonging 
to the secondary peak and lying between the two peaks re- 
spectively. Dark red, green and blue lines in the lower panel 
reproduce the locu s of red sequenc e , gree n valley and blue 
cloud as defined in ISalimbeni et all (I2OO8I ). 

We expect the two scenarios we discuss in the previous 
paragraph to predict a different colour evolution for galax- 
ies. In scenario A, we assume a continuity between the colour 
of galaxies belonging to the primary peak, valley and sec- 
ondary peak, in terms of blue cloud, "green valley" and red 
sequence. Since scenario A assumes a continuous motion of 
galaxies from the blue cloud to the red sequence, we expect 
galaxies with different SFRs to populate well defined regions 
of the colour-magnitude diagram. The relative abundance of 
galaxies in the different regions of the diagram is an indica- 
tion of the time they spend in the different phases. 

On the other hand, in scenario B we expect a larger 
spread of galaxies belonging to the secondary peak in the 
colour magnitude diagram. The physical processes responsi- 
ble for the retriggering have a strong dependence on the 
stellar mass and/or environment of galaxies, and are al- 
most insensitive to the position of the galaxy in the colour- 
magnitude diagram at the time of the retriggering. There- 
fore, the relative tightness of the locus where secondary peak 
objects are found in the rest-frame colour-magnitude dia- 
grams is in qualitative agreement with the former scenario. 
On the other hand, galaxies with SFR intermediate between 
the two peaks show a much wider distribution in fig. |4] We 
warn the reader that these are qualitative conclusions that 
rely on the assumption that the U — V colour is a fair tracer 
of the SFR activity, and that this colour is also sensitive 
enough to the low-levels of SFR for galaxies in the secondary 
peak. The fact that most of SFR estimates for sources below 
the SFRF peak have been derived from SED fitting also pre- 
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vents firm conclusions on this issue from GOODS-MUSIC 
data alone. A more quantitative analysis will require bet- 
ter estimates for SFR levels in galaxy population and better 
tracers of physical state of galaxies belonging to the different 
regions of the colour-magnitude diagram, and it is beyond 
the aims of the present work. 



5 CONCLUSIONS 

In this paper we consider the GOODS-MUSIC sample to 
study the SFR function and its redshift evolution in the 
redshift range 0.4 < z < 1.8, just after the peak in the cos- 
mic SFR, a fundamental epoch to understand the physical 
processes responsible for the decline of star forming activity 
as a function of galaxy stellar mass. In order to reduce the 
effect of completeness corrections, we restrict our analysis 
to galaxies with M* > 10^^ Mq, since at these stellar masses 
the GOODS-MUSIC is expected to be a complete sample. 
We thus compute the SFRF using SFR estimates coming 
from both SED fitting and/or a combination of UV light 
with 24/im observations. We compare the resulting SFRF 
with the predictions of 3 independent se mi-analytical mod- 
els of galaxy formati on and evolution l|Wang et all l2008l : 
ISomerville et al.|[2008l and morgana) and we obtain the fol- 
lowing results: 

• The agreement between the observed and predicted 
high-SFR end (SFR> lOMg/i/r) of the SFRF is good over 
the whole redshift range (once the observational errors are 
taken into account). This result implies that theoretical 
models are able to reproduce the space density evolution 
of the strongest starbursts. It is worth noting that the pre- 
dictions of the three models are remarkably similar in this 
SFR range. 

• The SFRF of massive galaxies, as seen by the GOODS- 
MUSIC is characterised by a complex evolution both in the 
typical SFR and in the normalisation. Theoretical models 
are able to cope with the SFR* evolution, but they predict 
a negligible density evolution. In particular, all models over- 
predict the space density of SFR ~ IMQ/yr galaxies. We in- 
terpret this behaviour as due to the well documented excess 
of inter mediate mass galaxies (lO^^M© < M* < 10^^ Mq) in 
SAMs (|Fontanot et al.ll2009bl '). 

• The observed low-SFR end of SFRF is characterised by 
a double peaked shape. Despite the uncertainties in the de- 
termination of such low SFRs from SED fitting methods, 
it is worth noting that none of the theoretical model we 
consider is able to reproduce this feature, and all of them 
predict a smooth decrease of the space density of low-SFR 
galaxies, with a well defined power-law slope a. Since this is 
also the SFR range where models differ most, each of them 
predicting a different value for a, we stress that future obser- 
vations providing stronger constraints on the low-SFR end 
of the SFRF will be of fundamental importance to under- 
stand the physical mechanisms responsible for the decline of 
SFR since z ~ 2 in massive galaxies. 

The main result of this paper lies in the analysis at the 
high-SFR end, since this region is completely sampled by 
24^m observations, which are critical for the recovery of the 
SFR levels. Also the evolution of the typical SFR* is well 
sampled by means of IR data out to z ^ 1.0, and allow us to 



put strong constraint on this relevant discrepancy seen be- 
tween data and models. Unfortunately, for lower SFR levels 
this information is not currently available, and we have to 
rely entirely on the results of SED fitting. Our tests suggest 
that the use of SED fitting results does not introduce any 
obvious distortion in the SFR distribution at intermediate 
SFR values, but we could not completely exclude the dou- 
ble peaked feature in the SFRF to be due to systematics in 
the SED fitting algorithm. From optical photometry alone is 
not possible to uniquely associate SFR ~ XM^Iyr galaxies 
(the region were models and data differ most and roughly 
corresponding to the region between the two peaks of the 
SFRF) with known populations in the colour-magnitude di- 
agram, i.e. the green valley or the blue cloud. The red U — V 
colours, with a relatively small scatter, of galaxies belong- 
ing to the putative secondary peak provide only a qualita- 
tive support to a scenario where star formation in massive 
galaxies decrease at a slower pace with respect to theoreti- 
cal expectations. If confirmed, this double peaked shape of 
the SFRF would thus represent a critical discrepancy be- 
tween models and observations, since all SAMs we consider 
predict a smooth, featureless, single slope power-law shape 
for the SFRF. Modifications in the assumed star formation 
law, taking into account the relation between star formation 
and the atomic and molecular gas content of galaxies may 
alleviate this tension, since the properties of model galaxies 
with low SFR leve ls are very sensitive to this modeling layer 
(|Lagos et al.ll201ll . see e.g.). 

On the other hand, thanks to the new facilities (e.g. 
Herschel), future observations in the IR region, able to char- 
acterise the properties of the IR peak, will finally determine 
the SFR levels of these objects and the relevance of this fea- 
ture. This is a fundamental task, since the low-SFR end of 
the SFRF is the region where SAM predictions differ most 
and therefore provides potentially strong constraints to the 
different approaches to star formation and stellar feedback. 
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